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Abstract The parametrically driven damped nonlinear Schrodinger equa- 
tion serves as an amplitude equation for a variety of resonantly forced 
oscillatory systems on the plane. In this note, we consider its nodal soliton 
solutions. We show that although the nodal solitons are stable against 
radially-symmetric perturbations for sufficiently large damping coeffi- 
cients, they are always unstable to azimuthal perturbations. The corre- 
sponding break-up scenarios are studied using direct numerical simula- 
tions. Typically, the nodal solutions break into symmetric "necklaces" of 
stable nodeless solitons. 



1. Two-dimensional localised oscillating structures, commonly referred to as 
oscillons, have been detected in experiments on vertically vibrated layers of gran- 
ular material pQ, Newtonian fluids and suspensions |2I3| . Numerical simulations 
established the existence of stable oscillons in a variety of pattern-forming sys- 
tems, including the Swift-Hohenbcrg and Ginsburg-Landau equations, period- 
doubling maps with continuous spatial coupling, scmicontinuum theories and 
hydrodynamic models |3I4| . These simulations provided a great deal of insight 
into the phenomenology of the oscillons; however, the mechanism by which they 
acquire or loose their stability remained poorly understood. 

In order to elucidate this mechanism, a simple model of a parametrically 
forced oscillatory medium was proposed recently The model comprises a 
two-dimensional lattice of diffusively coupled, vertically vibrated pendula. When 
driven at the frequency close to their double natural frequency, the pendula 
execute almost synchronous librations whose slowly varying amplitude satisfies 
the 2D parametrically driven, damped nonlinear Schrodinger (NLS) equation. 
The NLS equation was shown to support radially-symmetric, bell-shaped (i.e. 
nodeless) solitons which turned out to be stable for sufficiently large values 
of the damping coefficient. These stationary solitons of the amplitude equation 
correspond to the spatio-temporal envelopes of the oscillons in the original lattice 
system. By reducing the NLS to a finite-dimensional system in the vicinity of 
the soliton, its stabilisation mechanism (and hence, the oscillon's stabilisation 
mechanism) was clarified [5]. 

In the present note we consider a more general class of radially-symmetric 
solitons of the parametrically driven, damped NLS on the plane, namely soli- 
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Figurel. Solutions of eq.Q: TZq(t) (thin continuous line), Hi(r) (thick line), 
n 2 (r) (dashed). 

tons with nodes. We will demonstrate that these solitons are unstable against 
azimuthal modes, and analyse the evolution of this instability. 

2. The parametrically driven, damped NLS equation has the form: 

itpt + V 2 V + 2\ip\ 2 ip -^ = h^j* - i-ytp. (1) 

Here V 2 = d 2 /dx 2 + d 2 /dy 2 . Eq.([Q) serves as an amplitude equation for a wide 
range of nearly-conservative two-dimensional oscillatory systems under paramet- 
ric forcing. This equation was also used as a phenomenological model of nonlinear 
Faraday resonance in water [3] ■ The coefficient h > plays the role of the driver's 
strength and 7 > is the damping coefficient. 

We start with the discussion of its nodeless solitons and their stability. The 
exact (though not explicit) stationary radially-symmetric solution is given by 

V>o = Ae-' w TZ (Ar), (2) 

where r 2 — x 2 + y 2 , 

A 2 = 1 + ^h 2 - 7 2 , 6 = i arcsin (Tj , 

and lZo(r) is the bell-shaped nodeless solution of the equation 

K rr + ~K - K + 2TZ 3 = (3) 
r 

with the boundary conditions lZ r (0) — 7£(oo) = 0. (Below we simply write 1Z 
for Ho.) Solutions of eq.© are well documented in literature see Fig^ 

3. To examine the stability of the solution @ with nonzero h and 7, we 
linearise eq.Q in the small perturbation 



<fy(x, t) = e (^- r )*- ie ± [ u (x) + w( X )], 
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Figure2. Stability diagram for two-dimensional solitons. The (7, h — 7)-plane is 
used for visual clarity. No localised or periodic attractors exist for h < 7 (below the 
horisontal axis). The region of stability of the soliton i/jq lies to the right of the solid 
curve marked "n = 0". Also shown are the regions of stability of the solitons ipi 
and ip2 with respect to the radially-symmetric perturbations. (These lie to the right 
of the corresponding curves in the figure.) 



where x = Ax., t = A 2 t. This yields an eigenvalue problem 

L\u = -(/! + r)v, (L — e)v = {[i — r)u, 
where r = 7/ A 2 and the operators 



U = -V 2 



l-2K 2 {f), L x = L Q -4K 2 (?), 



(4) 



(5) 



with V 2 = d 2 /dx 2 + d 2 /dy 2 . (We are dropping the tildas below.) For further 
convenience, we introduce the positive quantity e = 2^/h 2 — j 2 /A 2 . Fixing e 
defines a curve on the (7, /i)-plane: 



Introducing 



h=^e 2 /(2-e) 2 + 1 2 . 

a 2 = M 2 - r 2 



(6) 
(7) 



and performing the transformation [§] 

«(x)^(/i + r)A- 1 V (x), 
reduces eq.Q) to a one-parameter eigenvalue problem: 
[Lq — e)v = \u, L\u = —Aw. 



(8) 



Wc first consider the stability with respect to radially symmetric perturba- 
tions u — u(r), v — v(r). In this case the operators (5J become 



dr 2 



1 d 

r dr 



+ l-2K 2 (r), Li = L Q -4K 2 {r). 



(9) 
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In the absence of the damping and driving, all localised initial conditions in the 
unperturbed 2D NLS equation are known to either disperse or blow-up in finite 
time |6I7I8| . It turned out, however, that the soliton i/jq stabilises as the damping 
7 is increased above a certain value [5]. The stability condition is 7 > 7 C , where 

2 ReA(e)ImA(e) 
7e = 7e(e) = — V(ImA)2 _ (R ^ A)2 - (10) 

We obtained A(e) by solving the eigenvalue problem JHJ directly. Expressing e 
via 7 C from (fTUfl and feeding into ©, we get the stability boundary on the 
(7, /i)-plane (Fig©. 

4. To study the stability to asymmetric perturbations we factorise, in (JSJ), 

u(x) = u{r)e mip , u(x) = v(r)e tmv , 

where tan<^ = y/x and m is an integer. The functions u{r) and v(r) satisfy the 
eigenproblem © where the operators © should be replaced by 

4 m) ^L + m 2 /r 2 , 4 m) ^i 1+ m 2 /r 2 , (11) 

respectively. This modified eigenvalue problem can be analysed in a similar way 
to eqs.©. It is not difficult to demonstrate that all discrete eigenvalues of (© 
(if any exist) have to be pure imaginary in this case, and hence the azimuthal 
perturbations do not lead to any instabilities of the solution in question jjjj. 

5. Besides the nodeless solution TZo(r), the "master" equation (© has so- 
lutions 1Z n (r) with n nodes, n > 1. (See Fig^}. These give rise to a sequence 
of nodal solutions ip n of the damped-driven NLS ©, defined by eq.© with 
TZq T^n- To examine the stability of the ip n , we solved the eigenvalue problem 
(jHl numerically, with operators as in (fTT|) . The radial stability properties 
of the nodal solitons turned out to be similar to those of the nodeless soliton 
ipQ. Namely, the ip n solutions are stable against radially-symmetric perturba- 
tions for sufficiently large 7. The corresponding stability regions for ipi and i\>i 
are depicted in Fig.©. However, the azimuthal stability properties of the nodal 
solitons have turned out to be quite different. 

Both ipi and ip2 solutions do have eigenvalues A with nonzero real parts for 
orbital numbers m > 1. (See Fig|3]and0]) Having found eigenvalues A for each 
e, one still has to identify those giving rise to the largest growth rates 

v = Rcfi-r (12) 

for each pair (e,7) [or, equivalently, for each (/i, 7)]. In (|12f) . /1 is reconstructed 
using eq.J7J|. The selection of real eigenvalues is straightforward; in this case we 
have the following two simple rules: 

• If, for some e, there are eigenvalues Ai > 0, A2 > such that Ai > A2, then 
v\ > V2 for this e and all 7. That is, of all real eigenvalues A one has to consider 
only the largest one. 

• If, for some e, there is a real eigenvalue Ai > and a complex eigenvalue A2, 
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Figure3. The discrete eigenvalues of the linearised operator JHJ for the one-node 
soliton, ipi. Panels (a) and (b) show the complex eigenvalues, Im A vs Re A. Arrows 
indicate the direction of increase of e. Panel (c) shows the real eigenvalues, as 
functions of e. 



with RcA2 > and Ai > RcA2, then v\ > V2 for this e and all 7. That is, one 
can ignore all complex eigenvalues with real parts smaller than a real eigenvalue 
— if there is one. 

The comparison of two complex eigenvalues is not so straightforward. In 
particular, the fact that ReAi > ReA2 does not necessarily imply that v\ > V2- 
Which of the two growth rates, v\ or is larger will depend on the imaginary 
parts of A, as well as on 7. 

In figures and we illustrate the real and imaginary parts of the eigen- 
values, arising for different m, for the solitons ip\ and ip%. The soliton ipi has 
discrete eigenvalues A associated with orbital numbers m — 0,1,..., 5 and the 
soliton tp2 with m = 0, 1, 10. 

In order to compare the conclusions based on the linearised analysis with 
direct numerical simulations of the unstable solitons ipi and ip2, we fix some h 
and 7 and identify the eigenvalue with the maximum growth rate in each case. 
In the case of the soliton we choose 7 = 3.5 and h = 3.6; the corresponding 




Figure4. The discrete eigenvalues of the linearised operator © for the two-node 
soliton, -02- Panels (a),(b), and (c) show the complex eigenvalues, ImA vs Re A. 
Panel (d) shows the real eigenvalues, as functions of e. 



e = 0.9146. The real and imaginary parts of A for each m as well as the resulting 
growth rates v are given in Table 1 (left panel). The eigenvalue with the largest 
ReA is associated with m = 0; however, for the given e and 7 the resulting 
v < 0. (This is because we have chosen a point in the "radially stable" part of 
the (7, /i)-plane, to the right of the "n = 1" curve in Fig|2) On the contrary, the 
growth rates corresponding to the real eigenvalues associated with to = 3,4, 5 
are positive for all 7. The maximum growth rate is associated with to = 4. 
The corresponding eigenfunctions u(r) and v(r) have a single maximum near 
the position of the minimum of the function lZ\(r); that is, the perturbation is 
concentrated near the circular "valley" in the relief of ip (x , y) | 2 . This observation 
suggests that for 7 = 3.4 and h = 3.5, the soliton ip\ should break into a 
symmetric pattern of 5 solitons ipo- one at the origin and four around it. 

Next, in the case of the soliton -02 we fix 7 = 4.5 and h = 4.53; this gives 
e = 0.6846. The corresponding eigenvalues, for each to, are presented in Table 1 
(right panel) . Again, the eigenvalue with the largest ReA is the one for to = but 
the resulting v is negative. The largest growth rates (1^3 = 0.54 and ^4 = 0.53, 
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Figure5. Evolution of the azimuthal instability of the one-node soliton. (a): the 
initial condition, soliton t/ji ; (b) and (c): dissociation of the ring-like "valley" into 
4 nodeless solitons; (d): divergence of the fragments. Here 7 = 3.5 and h = 3.6; 
shown is Reip. Note the change of the vertical scale in (d). 



respectively) are those pertaining tom = 3 and m = 4. The corresponding eigen- 
functions have their maxima near the position of the minimum of the function 
IZ2 (r). Therefore, the circular "valley" of the soliton ?/>2 is expected to break into 
three or four nodeless solitons ipo- (Since V3 is so close to 1/4, the actual number 
of resulting solitons — three or four — will be very sensitive to the choice of 
the initial perturbation.) Next, eigenfunctions pertaining to m = 5, 6, ...10 have 
their maxima near the second, lateral, maximum of the function TZ^if). The 
largest growth rate in this group of eigenvalues arises for m = 8. Hence, the 
circular "ridge" of the soliton ip2 should break into 8 nodeless solitons, with this 
process taking longer than the bunching of the "valley" into the "internal ring" 
of solitons. 

The direct numerical simulations corroborate the above scenarios. The tp\ 
soliton with 7 = 3.5 and h — 3.6 splits into a constellation of 5 nodeless solitons: 
one at the origin and four solitons of opposite polarity at the vertices of the square 
centered at the origin. The emerging nodeless solitons are stable but repelling 
each other, see Fig|3J Hence, no stationary nonsymmetric configurations are 
possible; the peripheral solitons escape to infinity. The tp2 soliton with 7 = 4.5 
and h — 4.53 has a more complicated evolution. As predicted by the linear 
stability analysis, it dissociates into 13 nodeless solitons: one at the origin, four 
solitons of opposite polarity forming a square around it and eight more solitons 
(of the same polarity as the one at the origin) forming an outer ring. (The fact 
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Figure6. The evolution of the azimuthal instability of the two-node soliton ip2- (a): 
the initial condition; (b)-(c): the rapid dissociation of the "valley" into 4 nodeless 
solitons; (c)-(d): a slower decay of the "ridge" into 8 solitons ip Q ; (e)-(f): the annihi- 
lation of the internal ring and the central soliton, and the repulsion of the persisting 
8 solitons. Here 7 = 4.5 and h — 4.53; shown is Reip. Note the change of the 
vertical scale in (e)-(f) w.r.t. that in (a)-(d). 
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Tablel. Eigenvalues A and corresponding growth rates v for the solitons tpi (left 
panel) and %l>2 (right panel). We included only the eigenvalues which can, poten- 
tially, give rise to the largest growth rate in each "symmetry class" m. Some other 
eigenvalues have been filtered out using the above selection rules. 



that the inner ring consists of four and not three solitons, is due to the square 
symmetry of our domain of integration which favours perturbations with m = 4 
over those with m = 3.) In the subsequent evolution the central soliton and the 
nearest four annihilate and only the eight outer solitons persist. They repel each 
other and eventually escape to infinity, see FigEI 

In conclusion, our analysis suggests the interpretation of the nodal solutions 
as degenerate, unstable coaxial complexes of the nodeless solitons ipo. 
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